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Hole-doped LaMnC>3 and electron-doped CaMnC>3 form self-trapped electronic states. The spectra 
of these states have been calculated using a two orbital (Mn e g Jahn- Teller) model, from which 
the non-adiabatic optical conductivity spectra are obtained. In both cases the optical spectrum 
contains weight in the gap region, whose observation will indicate the self-trapped nature of the 
carrier states. The predicted spectra are proportional to the concentration of the doped carriers in 
the dilute regime, with coefficients calculated with no further model parameters. 



I. INTRODUCTION 

Qualitative discussions of optical absorption from po- 
larons [jlj have often been given. For the heavily doped 
manganites such as Lao.7Cao.3Mn03, polarons form at 
temperatures T > T c ss 250K where ferromagnetic order 
is lost. Dramatic thermal shifts of optical conductivity 
agree j| with the qualitative polaron picture [Q. In 
the lightly doped end-member compounds LaMnOs and 
CaMnOa, polarons must exist at low temperature in the 
spin and orbitally ordered states. We have previously 
calculated the properties of these polarons in the 
Mn e g Jahn- Teller (MEJT) model. Here we point 
out that the polarons have a rich spectrum of local ex- 
cited states. We provide a quantitative calculation of the 
corresponding polaron-induced mid-infrared absorption. 
The transitions between these local levels are broadened 
by the Franck-Condon effect since each local level prefers 
its own rearrangement of the lattice. This is a general 
phenomenon expected in many insulators and may ac- 
count for mid-infrared absorption in other systems as well 

In the usual theory of small polaron absorption ||1|,|12| , 
an electron is trapped at a single site in a non-degenerate 
orbital. The photon permits a transfer of the electron 
to a first neighbor site. Direct transfer to the gound 
state first neighbor polaron (corresponding to dc con- 
ductivity) is exponentially suppressed by unfavorable vi- 
brational overlap ("Huang-Rhys") factors. The largest 
vibrational overlap is for transitions into states whose 
maximum amplitude overlaps with the zero-point motion 
of the distorted polaron on the original site. A Gaus- 
sian envoelope of Franck-Condon vibrational sidebands 
is predicted. The cases considered here differ in three 
ways. (1) Our electron levels are doubly degenerate be- 
fore distortion. (2) Hund's energy is large which blocks 
hopping unless the neighbor atom has spin ferromagnet- 
ically aligned. (3) Our polarons are not strictly confined 
to a single site, and consequently have localized excited 
states in the self-organized potential well. These com- 
plications are not all unique to manganites, and should 
affect polaronic absorption in other materials as well. 

Pure LaMnOs has a single Mn d electron in the e g dou- 



blet (V> 2 = {x 2 - y 2 )/V3, -03 = 3z 2 - r 2 ). Therefore there 
is a cooperative Jahn- Teller effect resulting in orbital or- 
dering. The peak |l3|jl4[| near 2 eV in <j(ll>) is interpreted 
]l4],p| as e g — > e g transitions across the Jahn- Teller (JT) 
gap. Pure CaMn03 has no e g states occupied, and a(u) 
peaks at 3.6 eV |Q due to excitations across the gap be- 
tween occupied (O 2p and Mn i2g|) states and empty (Mn 
e g |) states. These materials are insulating antiferromag- 
nets (AF) with the nice property of being continuously 
dopable. Here we focus on two cases, x = 1 — e (electron- 
doped CaMnOs) and x = e (hole-doped LaMnOs.) Both 
materials are insulators when pure, and remain so when 
lightly doped because of polaron formation. 

We use the model Hamiltonian Hmejt = Ht + Ti. C n + 
T-th + Hi/ § f§] similar to the one used by Millis [fLO), 
except that we have U = oo for LaMnOs. For CaMn03 
we include also the t2 g AF exchange term Tij and do 
not need the Hubbard U term for dilute doping. Two 
e g orbitals per Mn atom are kept to fully represent the 
symmetries of the orbitals and the crystal, and to drive 
JT distortions. In the hopping term 
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there is only one parameter, namely, the nearest-neighbor 
t = (dda) integral from Slater-Koster two-center theory. 
In the electron-phonon interaction, wc include both the 
oxygen breathing vibrations Q\ and the local -E s -type 
oxygen motions {Q2, Q3) which give the Jahn- Teller dis- 
tortion. 
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V^Ql(i)/a/s). (2) 



The indices a and (3 run though orbitals (^2,^3)- T~Cl is 
the lattice elastic energy of oxygen displacement along 
Mn-O-Mn bond directions with an Einstein restoring 
force. Note that the JT modes Q a (l) and Q a (V) are 
not independent when Mn sites I and I' are first neigh- 
bors, since the intervening oxygen is shared by both Mn 
atoms (this is called the "cooperative" phonon effect.) 
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At the LaMn03 end of the phase diagram, the domi- 
nant term of the Hamiltonian after electron-electron re- 
pulsion is H cp - Instead of delocalizing in a conduction 
band, holes in lightly doped LaMnC>3 lower their energy 
more by self-trapping with local oxygen rearrangement, 
forming "anti-JT" small polarons ||. At the CaMnC>3 
end, at the cost of AF exchange energy, spin-flipping 
in the AF background allows the doped-in electron to 
lower its energy by derealization. This spin-polaron ef- 
fect competes with a localizing JT polaron effect from 
Ti cp . The competition determines the size of the spin- 
lattice polaron jo]. 

In this work we calculate the excitation spectrum of the 
spin-lattice polaron in CaMnC>3 by Born-Oppenheimer 
approximation, and for LaMn03 we calculate the excita- 
tion spectrum of the self-trapped hole in a Cho-Toyozawa 
||l5f approximation. Our non-adiabatic solutions allow us 
to calculate the Franck-Condon broadened optical excita- 
tion spectra of these lightly doped materials and predict 
new spectral features that should be experimentally vis- 
ible. 



II. SPIN-LATTICE POLARON SPECTRUM IN 
DOPED CAMNO3 

Undoped CaMn0 3 has Mn 4 +(^ g ) ions with S = 3/2 
on an approximately simple cubic lattice, and bipartite 
(G type) AF order below Tn — 125K. One single doped- 
in electron will align with the S — 3/2 core spin by Hund 
coupling. Hopping to first neighbors is then blocked by 
the AF order. To reduce energy by derealization, one 
Mn S — 3/2 atom flips. Then the nearby oxygens distort, 
to lift the e g degeneracy. This gives a seven-site spin- 
lattice polaron Q. The parameters were chosen to be 
\t\/JS 2 = 158, and T = g 2 / K\t\ = 0.25. Each of the 
seven Mn atoms provides two e g orbitals to form a 14- 
dimensional Hilbert space in which the doped-in electron 
is partially delocalized. 

The ground state of this polaron is calculated in adi- 
abatic approximation. With the local oxygens distorted 
in the optimized pattern, the ground state lies in the gap 
with energy lowered from the e g band by — 2.13|i|. The 
remaining 13 excited states contain three dipole-allowcd 
optical excitations. By re-optimizing the oxygen distor- 
tion pattern, these states (with original energies 1.20|t| 
and 1.93|t| above the e g band) can have their energy low- 
ered to sit below the e g band by — 0.67|i|. When the elec- 
tron is optically excited to one of these states, oxygens do 
not have time to rearrange to a new optimal pattern, and 
therefore are left in a combination of vibrationally excited 
states with respect to the re-optimized oxygen pattern of 
the electronic excited state. The Franck-Condon prin- 
ciple therefore applies and gives vibrational sidebands. 
The energy level diagram is shown in Fig. |l|. 

To examine the optical spectrum of this seven-site po- 
laron surrounded by 36 oxygens, we start from one of 



the three symmetry-equivalent electronic ground states 
ipo(r,u ) (e.g., 8 = (|]) with local lattice distortion 
Uo = ( u o,i, u o,2, uo,36) which minimize the energy 
Eq(u) of the ground state. For distortions it near uq, the 
vibronic eigenstates associated with the electron ground 
state are approximated by the Born-Oppenheimer prod- 
uct: 

36 

*0,n(r,if) =^0(^)^0) Y\_4>i,m{Ui-Uo,i), (3) 
i=l 

where 4>i^ ni [ui — Uo,i) is the ni phonon state of the har- 
monic oscillator centered at Mo,i ■ The usual approxima- 
tion is made here that ipo(r,u) can be accurately rep- 
resented by ipo(r,uo). The true electronic ground state 
varies slowly with u, and also depends on the direction 
u — uq of deviation from ud, but corrections are small 
at low temperatures where the initial state is not vibra- 
tionally excited. Then the adiabatic ground state energy 
sheet E (u) depends quadratically on (ui — UQ, :i ) and we 
treat the oxygens as 36 independent Einstein oscillators, 
with energies 

36 

E , n = E (u ) + huj ^(n t + -). (4) 
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FIG. 1. Energy level diagrams for self-trapped states, rela- 
tive to the lowest energy of a state in an undistorted environ- 
ment. On the left are electron states in CaMnOa on the right 
are hole states in LaMn03. The lowest level is the ground 
state of the doped in e g electron, bound by spin and lattice 
distortions. The dipole-allowed excited levels shown vertically 
keep the ground state distortions (dashed lines for E\\z and 
solid lines for £ 1 z). Shown to the right are the relaxed 
energy levels corresponding to these excited states, and the 
shaded regions represent vibrational excitations accompany- 
ing the optical transitions. 
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Among the 13 polaronic excited states, there are only 
three states that couple to the ground state by light. 
They are \tp z >= (\ip3,z > —\ip 3 ,—z >)/V% an d two 
other states which are rotations of \ip z > to x and y di- 
rections. Here |^3, z > denotes the tp3 orbital of the Mn 
atom sitting immediately above the central Mn. These 
three excited states have their own optimal lattice distor- 
tion patterns u z , u x and u y . Starting from u~o, these elec- 
tronic states remain eigenstates along the paths towards 
u z , u x and u y in the lattice distortion space. Therefore 
in these chosen directions, the vibronic eigenstates asso- 
ciated with the electronic excited states have the product 
form, e.g.: 

36 

ty z ,n(r,u) = i)z(r,u z )Yl4>i tni (ui - u z> i). (5) 



Consider the optical excitation from \l/o to ^ z for ex- 
ample. The T = OK spectrum is given by 



(6) 



The matrix element has both an electronic part and a 
phonon part. The electronic part involves the momen- 
tum matrix element < ip z \p z 1 1/^3 , >, where \ip3,0 > de- 
notes the "03 orbital of the central Mn atom. This matrix 
element is < ip z \p z \4>3, >= ima/h < V'z|'^t| ? /'3i >j 
where m is the electron mass and a the lattice constant. 
This result can be derived by looking at FM spin order, 
and requiring that the corresponding FM Bloch eigen- 
states \k > of Ht (with energy e(fc)) should have velocity 
< k\p z /m\k >=de(k)/dhk z . 

The phonon part of the matrix element is |< 
nf, u z \n, uq >| 2 , where the initial lattice vibrational state 
has n vibrational quanta relative to the oxygen equilib- 
rium positions uq. The lattice vibrational state is un- 
changed during the electron excitation, but must be ex- 
panded in the new basis (denoted by quanta n') around 
the new oxygen equilibrium positions u z . We consider 
only zero temperature optical conductivity, which re- 
stricts the initial vibrational quantum numbers n to be 
zero. The non-zero overlaps of this initial phonon state 
with the new phonon basis states of the displaced lattice 
give a spectrum with a Gaussian profile of delta function 
peaks of the multi-phonon transitions accompanying the 
electron exitation. The phonon frequency is taken to be 
hui = 0.084 eV. This is a little higher than the stretching 
frequency reported from Raman data [fl6|| , but the dis- 
crepancy does not affect our conclusions. The hopping 
integral t is taken from density-functional calculations 
p7[ to be t = —0.75 eV. The delta functions are broad- 
ened into Gaussians of width \Zn/j where 7 = O.OMeV 
and n' is the total vibrational quantum. The results are 
shown in Fig. 0. 
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FIG. 2. Polaron-induced optical conductivity spectra of 
lightly doped CaMnOa. (a) The dashed curve is for E _L z 
and the solid line for E || z. Vertical lines present the adi- 
abatic solution under fixed oxygen distortions and show the 
peaks at 2.50 eV and 3.05 eV. (b) In (a), the solutions belong 
to one of the possible ground states (6* = 0) Q . There are two 
other possible ground states with orientations 9 = 27r/3, and 
6 = 47r/3, which permute the directions x, y and I. Therefore 
we show the average spectrum. The absolute value of cr(u}) is 
shown for 10% electron doping. 



III. POLARONIC SPECTRA IN HOLE DOPED 
LAMNO3 

In LaMnOs the electron-phonon term Ti. cp stabilizes an 
orbitally ordered ground state via a cooperative JT oxy- 
gen distortion. We use the same model as in CaMnC>3, 
with different filling, and fit it to the JT ground state. We 
previously used this model to predict the resonant multi- 
phonon Raman cross-section J^] and the angle resolved 
photoemission spectrum including the Franck-Condon 
broadening . The photoemission process creates a sud- 
den hole without time for oxygen relaxation. Thus it is 
in a vibrationally excited state. In the doped compound 
a hole is initially present in its ground vibrational state 
in a locally distorted oxygen environment. The zeroth 
order electronic many-body wave function has a hole h 
in the orbital at the origin (taken arbitrarily to be an 
x-oriented orbital) in an otherwise perfectly correlated 
state with one E g electron per atom. 



\h0> =c x (0)|JT> 



IJT > 



n4wii c ^')io> 



(7) 



and 



Here the orbitals tpx & n d i/V are orthogonal x- 
y— oriented orbitals (ip2 ± V^Vv^- 

At the next level of approximation the state \h0 > is 
coupled to 18 other states. Six of these states keep the 



3 



hole at the origin but have an orbital defect, or "orbiton" 
o on one of the surrounding Mn sites ±£ where £ is x, y, 
or z, 

\hO,o±i>=4 c (±i)c Y (±i)\hO> (8) 

Six more states have the hole moved to the neighboring 
Mn site, with no orbital defect, 

\h±i>= c Y (±i)4 c (0)\hO > . (9) 

The last six states have the hole on a neighbor and leave 
a misoriented \Y > orbital at the origin, 

\h±!,o0 >= c Y (±i)c Y (0)\hO > . (10) 

The Hamiltonian 7Ymejt becomes a 19 x 19 matrix in 
this truncated Hilbert space. 

In the adiabatic solution of the zcroth order hole 
ground state \h0 >, the neighboring oxygens in 
±x,±y,±z directions are distorted to positions given in 
Table [j], row 1 , which minimize the polaron energy || . 

With the oxygen pattern frozen in these positions, the 
diagonal element < h0\H cp \h0 > is set to zero. Non- 
zero Q3 distortions on the neighboring sites couple the 
state \h0 > to the six orbital excitations \h0, o±£ >. 
Row 2 of Table | gives the diagonal elements < hO, o ± 
£\H cp \hO, o ± £ > and row 3 gives the corresponding off- 
diagonal elements. The electron hopping term Ti t couples 
the state \h0 > to states with the hole moved to the one 
of the six neighboring sites, both without (Eq. |9|) and 
with (Eq. [h]) an orbital defect. The energies of these 
states are given in Table [| rows 4 and 5. 

Diagonalization of this 19 x 19 matrix gives the ground 
state polaron solution and the excited spectra. Among 
the 18 excited states, there are 9 odd states which are 
coupled by light to the even parity ground state with 
momentum matrix elements < ^Ipl/iO >= ima/h < 
^f\H t \hO >, whereas the other 9 even states are optically 
silent. No new parameter is needed for dipole matrix ele- 
ments to predict the absolute value of ct(cj). The spectra 
would consist of 9 delta functions, three for £Mield po- 
larized in z direction erj_(ci;) and six for o"|i(w) shown by 
vertical lines on Fig. (||, a). The parameters were cho- 
sen to be A = 0.95 eV, t = 0.5 eV, fiu = 0.075 eV, and 
Mn-Mn distance a = 4.0A. 

In order to investigate the lattice-coupled electronic ex- 
citation, the surrounding oxygens should be allowed to vi- 
brate. Since the frozen oxygen positions are not optimal 
for the excited hole states, Franck-Condon vibrational 
sidebands should appear in the optical spectrum. As- 
suming that the individual vibrational sidebands are not 
energy-resolved, there will be broadening by an amount 
roughly equal to the energy gain due to structural relax- 
ation. In order to quantitatively describe the process, a 
theory beyond the adiabatic approximation is needed. 

We use an approach similar to that applied to describe 
photoemission in the undoped material M. The method 



was proposed by Cho and Toyazawa |15j, who used an 
infinitely large truncated Hilbert space to diagonalize a 
ID polaronic Hamiltonian. For each of the states ap- 
pearing in Eq. @, Eq. (||) and Eq. (|ili|), we re-optimize 
the positions of the 6 oxygens surrounding the hole. The 
relaxed energies are listed in Table |, row 6, 7, and 8. 
Notice that for each of the states in Eq. (g), the energy 
gain from relaxation is A, similar to exciton creation in 
the parent compound ||. Each of these states, including 
I hO > , is now allowed to have arbitrary number of vibra- 
tional quanta of the surrounding oxygens around their re- 
optimized positions for that particular state. This adds 
additional vibrational multiples ntko to the relaxed en- 
ergy. In this chosen basis, the Tt cp + "Hi term of the 
Hamiltonian is diagonal. The complete Hamiltonian with 
the off-diagonal term 7i t can be diagonalized exactly. 
The hole ground state is not vibrationally excited at zero 
temperature, and the electronic part of this ground state 
has even symmetry which is coupled by light to odd ex- 
cited states. 

The Cho-Toyozawa calculation can be simplified by 
ignoring the sharing of the oxygens by its two nearest 
Mn neighbors. This partial removal of the "cooperative 
phonon" constraints was found in our photoemission cal- 
culation Q] to make only a small error. The energies 
of the odd parity states can be obtained as roots of the 
following analytic equation: 

G\{E)Gl{E) + G\{E)G\{E) = l 

G\{E) = t\ exp (-oj) (ID 

n— £ 

The roots E give energies of the excited states for light 
polarized in £ = x,y,z directions. The parameters a\, 
b\ and t\ in Eq. |ll| are given in table [j| The meaning 
of parameters a\ is the energy gain due to the structural 
rearrangement. The parameters b l a are energies of the 
electronic excited states in vibrational ground state with 
relaxed oxygen distortions. The hopping matrix elements 
t l a couple the ground state to the electronically excited 
states and they are proportional to the dipole matrix 
elements of the allowed transitions. The parameters A 
and t used in Eq. |ll| are all in units of %u> = 75 meV. 

The energy of the ground state is lowered by 0.39 
eV due to the hole derealization on the neighboring 
sites such that only 83 % of the polaron is present on 
the central site. This value is very close to the estima- 
tion by second order perturbation theory j^] , which gives 
-0.167i/r-0.488rt = 0.41 eV (r = A/(8i)). The excita- 
tion spectrum is measured from the renormalized ground 
state energy and shown on Fig. ||. It should be noted 
that the electronically excited states have their maximum 
weight on the neighboring sites and are not allowed to 
lower their energies by spreading on the next neighbor- 
ing sites due to the truncation of the Hilbert space. We 
have done an estimate of the finite Hilbert space effect 
by comparing the adiabatic spectra using 19 electronic 
states on 7 sites with those from the 79 states on 27 sites 
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including the next neighboring atoms. The spectrum in 
a larger basis has been shifted by approximately 0.3 eV 
to lower energies. The energy peak positions and rela- 
tive intensities in the adiabatic solution follow the trend 
of the non-adiabatic result (see Figj3la). We expect an 
overall shift of the spectrum in Fig. gto the lower energy 
by about 0.3 eV, such that the the spectrum sets in at 
about 0.7 eV. 
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FIG. 3. Polaron-induced optical conductivity spectra of 
lightly doped LaMnOa. (a) The dashed line is for E \\ z 
and solid line for E _L z. Vertical lines is an adiabatic solu- 
tion (frozen oxygen distortions). The non adiabatic spectra 
consist of three and six Franck-Condon broadened peaks for 
E || z and E _L z respectively, (b) Shows the average over 
polarizations. The overall spectral weight is proportional to 
the hole concentration and is given for 10% doping. 



IV. CONCLUSION 

The particular details for the manganites show that 
interesting processes occur which were not fully antici- 
pated in the qualitative discussions previously given |Q. 
Similar processes should be at work in other insulators 
such as cuprates and nickelates. 

New spectral features are predicted to appear in the 
JT gap of lightly doped LaMnC>3, and in the insulating 
gap of lightly doped CaMnC>3. The excitation spectra are 
calculated in non-adiabatic Born-Oppenheimer or Cho- 
Toyozawa approximation. These excitations are dipole- 
allowed with dipole matrix elements proportional to the 
hopping parameter t. 

In the high T paramagnetic phase, in a mean field ap- 
proximation, spin disorder has the effect of reducing the 
parameter t by 30%. This estimate suggests that the 
spectral weight due to hole excitations in LaMnC>3, for 
example, is reduced by 50% above the Neel temperature 
Tjy ~ 140-ftT |18|, whereas the orbiton excitation cen- 



tered at 2 eV H is not sensitive to magnetic order and 
hence is temperature independent. In CaMnOa, above 
Tjm « 125K spin polarons should persist in disor- 
dered and increasingly delocalized form. This will cause 
the spectrum of Fig. || to broaden and shift to lower 
energy. 

Our picture makes the prediction that, at frequencies 
well below the lowest weak peak in the sub-gap region, 
there should be a weak Urbach tail with characteristic 
Urbach temperature dependence. The tail is from the 
lowest-energy processes, and are weak because of small 
vibrational overlap. The available experimental data 
on doped compounds agree with our picture. Further 
experiments on more lightly-doped materials, and espe- 
cially transmission experiments through thin single crys- 
tals, would give a better test of our theory. 



ACKNOWLEDGMENTS 

This work was supported in part by NSF Grant 
No. DMR-0089492 and by DOE Grant No. DE-AC-02- 
98CH10886. 



[1] D. Emin, Phys. Rev. B 48, 13691 (1993). 

[2] Y. Okimoto and Y. Tokura, J. Supercond. 13, 271 (2000). 

[3] K. H. Kim, J. H. Jung, and T. W. Noh, Phys. Rev. Lett 
81, 1517 (1998); S. G. Kaplan, M. Quijada, H. D. Drew, 
D. B. Tanner, G. C. Xiong, R. Ramesh, C. Kwon, and 
T. Venkatesan, Phys. Rev. Lett. 77, 2081 (1996); H. L. 
Liu, S. L. Cooper, and S.-W. Cheong, Phys. Rev. Lett. 
77, 4684 (1998). 

[4] A. J. Millis, P. B. Littlewood, and B. I. Shraiman, Phys. 
Rev. Lett. 74, 5144 (1995); A. J. Millis, R. Mueller, and 
B. I. Shraiman, Phys. Rev. B 54, 5405 (1996); S. Yunoki, 
A. Moreo, and E. Dagotto, Phys. Rev. Lett. 81, 5612 

(1998) ; A. S. Alexandrov and A. M. Bratkovsky, Phys. 
Rev. B 60, 6215 (1999). 

[5] P. B. Allen and V. Perebeinos, Phys. Rev. B 60, 10747 

(1999) . 

[6] Y. R. Chen and P. B. Allen, Phys. Rev. B 64, 064401 
(2001). 

[7] V. Perebeinos and P. B. Allen, Phys. Rev. Lett. 85, 5178 

(2000) . 

[8] P. B. Allen and V. Perebeinos, Phys. Rev. Lett. 83, 4828 
(1999). 

[9] V. Perebeinos and P. B. Allen, Phys. Rev. B 64, 085118 

(2001) . 

[10] A. J. Millis, Phys. Rev. B 53, 8434 (1996). 

[11] J. D. Perkins, J. M. Graybeal, M. A. Kastner, R. J. Bir- 
geneau, J. P. Falck, and M. Greven, Phys. Rev. Lett. 71, 
1621 (1993); D. A. Crandles, T. Timusk, J. D. Garret, 
and J. E. Greedan, Physica C 216, 94 (1993). 



5 



[12] A pedagogical discussion can be found at 
http:/ /insti. physics. sunysb.edu~allen/Polarons/MIR. html 

[13] T. Arima, Y. Tokura, and J. B. Torrance, Phys. Rev. B 
48, 17006 (1993). 

[14] J. H. Jung, K. H. Kim, D. J. Eom, T. W. Noh, E. J. 
Choi, J. Yu, Y. S. Kwon, and Y. Chung, Phys. Rev. B 
55, 15489 (1997); J. H. Jung, K. H. Kim, T. W. Noh, E. 
J. Choi, and J. Yu, Phys. Rev. B 57, 11043 (1998). 

[15] K. Cho and Y. Toyozawa, J. Phys. Soc. Jpn. 30, 1555 
(1971). 

[16] E. Liarokapis, Th. Leventouri, D. Lampakis, D. Palles, 
J. J. Neumeier, and D. H. Goodwin, Phys. Rev. B 60, 
12758 (1999); E. Granado, N. O. Moreno, H. Martinho, 
A. Garcia, J. A. Sanjurjo, I. Torriani, C. Rettori, J. J. 
Neumeier, and S. B. Oseroff, Phys. Rev. Letters 86, 5385 
(2001). 

[17] W. E. Pickett and D. J. Singh, Phys. Rev. B 53, 1146 
(1996). 

[18] E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 



TABLE I. The parameters of the single-hole version of 
Hamiltonian T^mejt for oxygen displacements (row 1), the 
diagonal elements < hO, o ± £\H\h0, o ± £ > are shown in row 
2, the off-diagonal coupling elements gQz(€) are in row 3, the 
diagonal elements < h ± £\H\h ± £ > are in row 4, the diago- 
nal elements < h ± I, oO > \H \h ± £, oO > are in row 5. The 
relaxed energies of each of the states |/i0, o ± I > are in row 
6, \h ± £ > are in row 7 and \h ± £, oO > are in row 8. The 
parameters a\, b\ and t\ (i = 1, 2, 3) used in Eq. [ll] are given 
in the rest of the table. 
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